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Abstract 

The formation of marginally trapped surfaces in the off-center collision of two shock waves on 
AdSx) (with D = 4, 5,6,7 and 8) is studied numerically. We focus on the case when the two 
waves collide with nonvanishing impact parameter while the sources are located at the same 
value of the holographic coordinate. In all cases a critical value of the impact parameter is 
found above which no trapped surface is formed. The numerical results show the existence of 
a simple scaling relation between the critical impact parameter and the energy of the colliding 
waves. Using the isometries of AdS^ we relate the solutions obtained to the ones describing the 
collision of two waves with a purely holographic impact parameter. This provides a gravitational 
dual for the head-on collision of two lumps of energy of unequal size. 
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1 Introduction 



The applications of AdS/CFT duality to the description of heavy ion collisions has attracted 
quite a lot of attention in recent years, in the hope of describing the phenomenology observed in 
RHIC and expected in the LHC (see pQ for reviews). The detection in RHIC of a large elliptic 
flow and jet quenching are interpreted as indications that the quark-gluon plasma produced in 
the collision is strongly coupled, thus making AdS/CFT a very useful tool to capture some of 
the relevant physics involved. In spite of the difference between real QCD and M = 4 SYM, 
the calculations of some relevant observables using gravitational duals |2j show a reasonable 
good agreement with the experimental results. 

A possible way to model the collision of two heavy ions in a strongly coupled gauge theory is 
to consider the scattering of two energy lumps in a M = 4 SYM theory at strong coupling j3l 0] . 
The gravitational dual of such a system is a space-time representing the collision of two shock 
gravitational waves propagating in AdSs. Nonetheless, the picture of the heavy ions proposed 
in [4] is not free of problems f5]. For example, it produces an energy density profile for the lumps 
in the gauge theory that falls off as a power of the transverse distance from its center, quite 
unlike the exponential fall off expected from phenomenological nuclear potentials. However, in 
spite of these caveats, these systems are a good test bench to analyze collective properties of 
strongly coupled plasmas that might be of relevance in understanding the phenomenology of 
heavy ion collisions [5]. 

In this paper we study various aspects of the collision of shock gravitational waves in AdS 
space-time in various dimensions. The metric in the region before the collision takes place (i.e., 
where at least one of the two null coordinates u or v is negative as shown in Fig. fll) is given by 



where $±(z,x) are the wave profiles and L is the AdS^ radius. A full description of the final 
fate of the wave collision (and its holographic interpretation) would require knowing in detail 
the solution to the Einstein equations in the interaction region u > 0, v > 0, something that up 
to date has not been fully achieved (for some important progress in this direction see [7]). Based 
on the intuition in flat space, however, one expects that the collision of two shock waves would 
produce some kind of dressed singularity. The holographic interpretation of such a process 
would be the thermalization of the Af = 4 SYM plasma produced as the result of the collision 
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of the two incoming lumps whose energy distribution is given by a CFT energy momentum 
tensor with components [8] 

W.W = vd(s D %p_ 6 > + ^] D -/ (t,) - (L2) 



with (;z±,&±) the holographic and transverse coordinates of the sources of the waves (1.1). 

A way to avoid solving for the geometry in the u > 0, v > wedge is to look for the 
formation of an apparent horizon before actually getting into the region of space-time where 
the interaction takes place [TTJ. In AdS for the case of zero impact parameter it was found in 
[I] that such a trapped surface is always formed for any energy of the incoming waves. The case 
of collisions with a nonvanishing impact parameter along the field theory coordinates has been 
studied in AdSs [TU] , where there is an energy-dependent critical value of the impact parameter 
above which no marginally trapped surface of the type sought is formed. The results found in 
[I] and [10] are qualitatively similar to the corresponding situations in flat spacetime [T2| [13]. 

In the dual Af = 4 SYM theory the interpretation of the results of [10J would be the existence 
of a critical impact parameter for the thermalization of the plasma following the collision of two 
energy lumps of the same size. This can be understood if we think that for large enough impact 
parameter the two energy distributions do not have enough overlap to induce a thermalization 
of the whole plasma once the collision takes places. Once again, to see what happens as the 
result of the collision would require solving the field equations into the interaction region^ 

Another interesting physical situation is that of the collision of two shock waves whose 



sources have different values of the holographic coordinates. A look at Eq. (1.2) shows that in 
the holographic theory this describes the collision of two energy lumps of different size, with or 
without impact parameter. This problem has been addressed in [8], where the trapped surface 
equation is solved analytically in the limit in which the impact parameter is much smaller than 
the critical value. It would be interesting, nevertheless, to study the collision of waves with 
"large" holographic impact parameter in various dimensions. 

In Ref. [H] the formation of trapped surfaces in the head-on collision of two shock waves 
was studied for incoming waves characterized by a finite size in transverse space, in contrast 



1 We should not forget that the non-existence of the trapped surface of Penrose's type [TT] does not exclude 
that other trapped surfaces are formed in the interaction region. 
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with the case studied in [4] where the energy density along the wave-front follows a delta- 
function distribution in transverse coordinates. For D = 4 and D = 5 it was found that the 
marginally closed trapped surface only forms if the transverse size of the wave is smaller than 
certain critical value. If we interpret the formation of the trapped surface as a signal of eventual 
black hole formation as a result of the collision, this could imply at face value the existence of 
a threshold for thermalization in the holographic gauge theory as a function of the spread of 
the wave in the gravitational dual. It is not clear, however, what meaning the spread of the 
gravitational sources introduced in [H] might have in the holographic conformal field theory. 
For example, a calculation of the holographic energy-momentum tensor shows that the size of 
the energy distribution in the boundary theory is independent of the value of the deformation 
parameter. Indeed, applying the holographic prescription [9] to the solution found in [H] for 
the wave profile x), the associated holographic energy-momentum tensor is 

2 D ~ 2 fi L D_1 

< T ->CFT = Vol (5^3) (f2 + jL2) D-2 5 W> ( L3 ) 

and similarly for (T w )cft by replacing u — > v . Thus no trace of the smearing of the gravita- 
tional source is left in the energy distribution in the boundary theory, its transverse size being 
determined solely by the value of the holographic coordinate of the source, z = L. 

In this note our aim is twofold: first to study numerically the formation of marginally closed 
trapped surfaces in the collision of two shock waves in AdS^ space-time with a nonvanishing 
impact parameter in the field theory coordinates, thus extending the analysis of [TU] to D ^ 5. 
Second, to apply the results obtained to the analysis of collisions of two AdS^ shock waves with 
"purely holographic" impact parameter, i.e. when the sources of the two incoming waves are 
located at points in transverse space that only differ in the value of the holographic coordinate. 
The strategy used consists in exploiting the underlying 0(2, D — 1) isometries of AdS^ to map 
this problem to a problem of the scattering of two waves with "spatial" impact parameter (cf. 
[SJ [TO])- As explained above, in the spirit of the AdS/CFT correspondence, collisions with 
pure holographic impact parameter provide the gravitational dual of the high-energy head-on 
collision of two energy distributions of unequal size in the strongly coupled CFT. 

In Section [2] we review the series of isometries of AdS^ that allow to rotate a general 
collision with purely holographic impact parameter into a symmetric collision with an impact 
parameter along the field theory coordinates. The numerical analysis of the trapped surface 
equations is done in Section [3] Finally, our results together with their holographic interpretation 
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Figure 1: The four regions in the space-time describing the interaction of two shock waves in 
AcISd- The blue (right) and red (left) thick lines represent the location of the incoming waves, 
supported respectively at u = 0, v < and u < 0, v — 0. 



been deferred to the Appendix. 

2 Gravitational shock wave collisions with holographic 
impact parameter 

We study the collision of two shock waves propagating on AdS^ whose delta-function sources 
are located in the hyperbolic transverse space Hd_ 2 at the points {z±,b±). This space can be 
seen as a two-sheeted hyperboloid 



in (D — l)-dimensional Minkowski space. Poincare coordinates (z, x) can be introduced 



are discussed in Section |4| To make the presentation self-contained some technical details have 
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in which the metric reads 

L 2 



ds* 



D-2 z 2 



(dz 2 + dx 2 y (2.3) 



The metric of the full space-time previous to the collision is of the form shown in Eq. (1.1 ) 
and is characterized by the profile functions <f>±(z,x). The Einstein equations reduce, outside 
the interaction region in Fig. [TJ to the Poisson-like equation in transverse space 

(□ Hd _ 2 - ^jj^j ^±{z,x) = -16nG N fi ± 5(z - z ± )5^\x -b ± ). (2.4) 

Here □h_ d _ 2 is the Laplacian in H jD _ 2 



r D-2 2 

- ' .4-Oa \ , V72 



□h d _ 2 = -j^-d, [z^d z j + -V 2 (2.5) 

and fi± is the energy of the incoming waves, as read from the energy-momentum tensor 

T uu = fi + ^Y' 2 6(u)5(z~z + )6^\x-b + ), 

T vv = »-(y) D ~ 2 S(v)5(z-z^ d ~ 3 \x-L). (2.6) 

In what follows we focus on the family of incoming trajectories considered in [3] where the 
sources are located at the points (z±, b±, 0, . . . , 0) with 

L 



z±(a) 



^l + /3 2 ± cos 



a 



, . . LB sin a , 

a/ 1 + a 2 ± p cos a 

with < a < 2n. 

The space Hd_2 has an 0(1, D — 2) group of isometries. This can be seen as the part of the 
full 0(2,D — 1) isometry of AdSu that preserves the null coordinates u, v. The group 0(1, D — 2) 
is generated by the 0(D — 2) rotations of the ambient space coordinates (Y 1 , . . . , Y D ~ 2 ) among 
themselves, plus boost along each of these "spatial" coordinates. In the following we describe 
a series of isometries of AdS^ (and therefore of M.0-2) that connect a general collision of two 
waves with sources at ( . L ,0,0,...) to a symmetric collision where the sources are located 
at (L, ±LB, 0, . . .). 
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0(2) rotations. To begin with we are going to focus on the 0(2) subgroup of 0(D — 2) that 
rotates the embedding coordinates Y 1 , Y D ~ 2 



a 



Y 

yiB-1 



Y 1 cos9 + Y D - 2 sin9, 
-Y 1 sinO + Y D - 2 cos6. 



(2.1 



In terms of the Poincare coordinates (z,x) = (z,x,Xt) the 0(2) rotations take the form 

2L 2 z 



x 



Si 



L 2 + z 2 + x 2 + x*t + (L 2 — z 2 — x 2 — Xj) cos 9 — 2Lx sin 9 ' 
2L 2 x cos 9 + L(L 2 - z 2 - xl) sin 6 



L 2 + z 2 + x 2 + x 2 + (L 2 



x 



2 — Xj) cos 9 — 2Lx sin 9 ' 



(2.9) 



2L 2 x T 



L 2 + z 2 + x 2 + Xj, + (L 2 — z 2 — x 2 — Xj) cos 9 — 2Lx sin 9 
In spite of the nonlinear action on the coordinates it can be checked without much effort that the 
previous transformations form an abelian group. Actually, the 0(2) rotations of the coordinates 
of M.o-2 can be written as the one-parameter group of transformations 



U{9) = exp 



(2.10) 



where 



-idj 

i ( z 2 + x 2 + x'i ) dj — 2ixj ( zd z + xd x + xt ■ V 



(2.11) 



are respectively the generators of translations and "special conformal transformations" on the 
field theory coordinates (x, xt)- This transformation can be lifted to the corresponding isometry 
of the full AdS/) space-time. 



Acting on the family of points shown in Eq. (2.7) the 0(2) rotations amount to a translation 
in the parameter a 

z' ± (9) = z± (a + 9), b' ± (9) = b ± (a + 9). (2.12) 

Hence, by performing an isometry of the background M.d-2 we can transform the points 

L 



z± 



0. 



(2.13) 



separated in the holographic direction z, into 

z' ± {9) = 

b' ± (6) = 



1 + /3 2 ±(3cos9' 
L/3 sm6 



(2.14) 



+ (3 2 ± (3 cos9' 

This 0(2) transformation maps a collision of two shock waves with holographic impact parame- 
ter into a problem in which the two waves collide generically with both spatial and holographic 
impact parameter. 



In our original problem (2.4) the group of isommetries 0(1, D — 2), and therefore its 0(2) 



subgroup (2.9), is broken only by the presence of the sources. Because of this one can use this 
underlying symmetry as a solution generating technique, since the transformation of a solution 
to the Einstein equations by an isometry is again a solution where the position of the sources 



is also tranformed. To see this we first compute the Jacobian associated with (2.9) 
d(z',x',x T ') T 2L 2 



d(z, x, xt) 



L 2 + z 2 + x 2 + x 2 + (L 2 - z 2 

D-2 



f? ) cos 9 — 2LX 1 sin 9 



D-2 



(2.15) 



Now we apply (2.9) to the equations (2.4) with sources located at (2.13). The Laplacian is 
invariant, □e c _ 2 = , since we are dealing with an isometry of Hd_ 2 . On the other hand, 

the transformation of the delta functions can be computed using the expression of the Jacobian, 
leading to the result 



$+(z',x") 



L 2 

\Q-kGnIj!±(9) 



4(g) 

L 



D-l 



6(z' - z' ± {9))s(x' - b' ± {9)y {D - A \x^). (2.16) 



This describes the collision of two AdS^ waves with sources at (|2.14|) and energies 



'l + /3 2 ±/3cos^ 



(2.17) 



^/T+P 2 ±f3 

This transformation of the energies fi± as the result of this rotation can be understood on 
physical grounds remembering that they correspond to the energy of the incoming field theory 
energy lumps at the boundary. The 0(2) transformation on the transverse hyperbolic space 
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brings one of the sources closer to the AdS boundary whereas the other "moves" deeper 
into the bulk (see Fig. [2]). From the point of view of the CFT on the boundary this means that 
the energy of the configuration associated with the first source increases whereas the energy of 
the lump associated with the second source gets smaller. This is a simple consequence of the 
IR/UV connection in AdS/CFT. 



We now particularize Eq. (2.16) to the case 9 = § , where the z coordinate of the two sources 



are equal. We find that the Einstein equations read in this case 



/ / \ D-l 

-167TGV4 f 5{z'- z' ± )d(x' - b' ± y^ D - 4 \x^) (2.18) 
and describes now the problem of the collision of two waves with energies 

/4 = ^ ± -Vl + £!_ (2.19) 

v/TT^±/3 

whose sources are located at 



z+ 



r L=, b' ± = ±^£^. (2.20) 

This means that we have transformed our problem into the collision of two waves with spatial 
impact parameter. 

Longitudinal boosts. This is however not the whole story. We have started with a collision 
where the two waves had different energy and, in general, this is also the case after the 0(2) 
rotation has been performed. For practical purposes, however, it is simpler to study symmetric 
collisions where the two waves approach each other at the same energy. To achieve this we can 
use another isometry of AdS/) consisting in longitudinal boosts on the u, v coordinates. To see 
how this affects the shock wave we apply the boost 



u 



Am, v = \~ 1 v, A>0 (2.21) 



to the shock wave metric (1.1). The metric preserves the same form now with the profile 

§±(z',x') = \ Tl <S>±(z',x'). (2.22) 
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Figure 2: Schematic representation of the action of the 0(2) rotation on the sources of the two 
colliding waves. 



Applying now the boost to the components of the energy-momentum tensor (2.6) we find that 
the energy n± associated with the transformed wave is 



— / \ =fi / 



(2.23) 



Thus, by applying to Eq. (2.18) a longitudinal boost with parameter 

we end up with a symmetric collision with energy y/JT^JH 
D — 2\ 



(2.24) 



/ / \ D-l 

16nG N ^jn ( Z f) 8(z' - z' ± js(x' - 6' ± )^- 4) (f T ')- 



(2.25) 



Coordinate rescalings. There is however something unpleasant about the expressions of z'^ 



± 



and b' ± shown in Eq. (2.20). Since our only free parameter is (3, we notice that it is not possible 
to change the spatial impact parameter b'± without at the same time moving the sources in 
the holographic direction z±. To consider collisions where the spatial impact parameter is 
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decoupled from the value of the holographic coordinate of the sources we apply the following 
rescaling to the coordinates 



z " = y /l + /3V, x" = y/1 + /3 V, x£ = y/T+Pxi. (2.26) 
In the full AdS^ space this coordinate transformation is accompanied by similar rescaling of 



the two null coordinates. In the new coordinates the sources (2.20) are located in the points 

(2.27) 



zl=L, 



thus decoupling the value of the holographic coordinate from that of the impact parameter. 
We have to keep in mind that this rescaling is an isometry of Hd_2. Implementing it in Eq. 



(2.25) we arrive at 



16nG N J^-5(z" - L)5(x" ± Lp)^ D - A \x^'). (2.28) 



1 + /3 2 

Now the two shock waves have spatial impact parameter b = \b'+ — b'!_\ = 2L/3 while both sources 
lie at z" = L. 

To summarize, we started with the collision of two shock waves with energies fi± whose 



sources were located at the points (z±,x±) given in Eq. (2.20). Taking advantage of the 
isometries of AdS^ (and of its transverse section Md-z) this problem has been mapped into the 



collision described by Eq. (2.28) where the two sources are located at z = L with an impact 
parameter b given by 



b = \b'L 



2L(3 



z\ = Az, 



(2.29) 



with Az the holographic impact parameter of the original collision. The change of coordinates 



connecting Eqs. (2.18) and (2.28) is given by 



L 2 + z 2 + x 2 + x*t + (L 2 — z 2 — x 2 — x*t) cos 9 — 2Lx sin 9 ' 



x 



2LV 1 + P 2 xcos9 + Ly/1 + f3 2 (L 2 - z 2 - x 2 ± ) sin^ 
L 2 + z 2 + x 2 + x| + [L 2 — z 2 — x 2 — x*t) cos 9 — 2Lx sin 9 ' 

2L 2 ^/l + f3 2 x T 
L 2 + z 2 + x 2 + x^ + (L 2 — z 2 — x 2 — Xj) cos 9 — 2Lx sin 9 
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(2.30) 



To find the relation between the wave profiles we recall that these functions are only transformed 
by the longitudinal boostJ^J Using the value (2.24) for the parameter of this boost we find that 



the solution $± (z, x) for the collision with holographic impact parameter can be written in 



terms of the solution Q±(z" , x") to Eq. (2.28) as 

$ ± (z,x) = ^ + * ± (Az,x),x"(z,x)), (2.31) 

^+ V v 1 + P - P I 



where the coordinate transformations z"(z,x) and x"(z,x) are given by (2.30). 

The bottom line of this whole discussion is that the solution to the problem of finding a closed 
trapped surface in the collision of two shock waves with an impact parameter aligned along 
the holographic direction is reduced to the one of finding this trapped surface in the collision 
of two shock waves with an impact parameter that is purely spatial. This latter problem is 
easier to tackle numerically. Once the trapped surface is found in this case the original one is 



obtained by applying the transformation (2.30). Because of the family of incoming trajectories 



(2.7) used our discussion applies to any collision problem where the position of the sources in 
the holographic coordinate z± satisfies z + Z- = L 2 . Notice, however, that for any other value 
of the product z + z_ a rescaling of all the coordinates in AdS can be performed to map it into 
the problem at hand. 

3 Numerical analysis of the formation of marginally closed 
trapped surfaces 

We have seen how the isometries of AdSi? can be exploited to rotate a collision of two shock 
waves with holographic impact parameter into a collision of two waves sourced at z = L and 
with spatial impact parameter. In D = 5 this latter problem has been solved numerically in 
[10J. In this section we extend the numerical solution to other values of the dimension and 
study the emergence of thresholds for the formation of the closed trapped surface as a function 
of the impact parameter. 

Before embarking on the search for closed trapped surfaces we have to switch to a new 
system of coordinates (U, V, Z, X) in which the null geodesies are continuous across the location 
of the waves [U [TU] (see the Appendix for the technical details). Once this is done, we look for 
2 We have to keep in mind that in AdS^ the rescaling ( 2.26| ) also affects the null coordinates. 
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<s 4 



(3-1) 



marginally closed trapped surfaces S that lie along the position of the colliding shock waves 
{U = 0, V < 0} U {U < 0, = 0} [TT1 H]. This surface can be written as the union of two 
branches S = S + U 5_ parametrized as 

u = -Mz,x) _f u = 

V = " W= -il>-(Z,jC] 

Although most of this surface lies outside of the region U > 0, V > where the two waves 
interact, it "feels" the interaction taking place through the matching conditions at the (D — 
3)-dimensional intersection of the two branches C = S + fl 5_ C Hd_ 2 . Imposing that the 
congruence of outer null geodesies normal to S has zero divergence leads to the equation (see 
[1] for the details) 

D-2' 



□if 



L 2 



± 



0. 



where Xj^) are the profiles of the incoming waves and 

* ± (z,x) = ^ ± (z,x). 



(3.2) 



(3-3) 



In addition, from the very form (3.1) of the sought marginally closed trapped surface and the 



requirement that the congruence of null geodesies is continuous across the intersection C we 
find the conditions 



*±(Z,X' 



with g a b the metric of 



0. 



g ab d a y ± (Z,X)d b V ± (Z,X] 



(3.4) 



In order to solve the problem (3.2), (3.4) we solve the differential equation numerically 



using the technique employed in Ref. [12] for the collision of shock waves in flat space-time. 
To implement the procedure we switch to radial coordinates in Md_2 defined in terms of the 
hyperboloid coordinates by 



Y° = r 2 + L\ Y^rd 1 



Y 



D-2 



D-2 



(3.5) 



where . . . , d D ~ 2 ) parametrize the unit sphere S D ~ 3 , i.e., (t? 1 ) + . . . + (i9 D ~ 2 ) 2 = 1. Now the 
metric takes the form 

dr 2 



+ r 2 d6 2 + r 2 sin 2 6dQ 



D-A- 



(3.6) 



L 2 
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Here we have singled out a polar angle 9 G [0, tt] on the sphere S D ~ 3 of constant r. In these 
new coordinates^] (r, 9, with i — 1, . . . , D — 3, the sources of the incoming waves are located 
respectively at the points (|, 0, 0, . . . , 0) and (|, 7r, 0, . . . , 0). Hence, the residual 0(-D — 3) 
symmetry is generated by the isometries of the sphere S D_4 at constant r and constant 9. 
As a consequence of this symmetry, both the profiles of the incoming waves &±(r, 9) and the 
functions defining the two branches of the trapped surface ^±(r, 9) only depend on r and 9. In 
addition, the surface C is generated by the action of 0(D — 3) on a curve that we parametrize as 
r = LG{9) for < 9 < n. We have introduced the AdS radius L in order to make the function 
G{9) dimensionless. 

Introducing the function H± (r 

2 \ _ fn-.-m.2_ _ , 

(3.7) 



1 + 



$±(r,i 

,2 



^±( r , #), Eq. (3.2) reads in these coordinates 
D - 



L 2 ) r rL 2 
with the boundary condition 



r 2 tan# 



da 



D-2 



H. 



± 







H ± (r,9) =$ ± (LG(6),0). (3.8) 

r=LG{6) 

This problem has to be solved in the region < r < LG(9), < 9 < tt, i.e. those points in the 
upper half-plane limited by the curve r = LG(9). Because the sources of the waves are located 
symmetrically around r = the function G(9) has the symmetry 



G(9) = G(ir 



whereas \l/±(r, 9) and $±(r, 9) satisfy 
\l/ ± (r, 9) = \I/ ± (r,7r - 



$ T (r, tc 



(3.9) 



(3.10) 



The boundary G(9) is determined by imposing that the solution to the boundary problem 
satisfies the addition constraint 

^,2 



1 + 



L 2 



1 



4. 



(3.11) 



r=LG{6) 

In the symmetric collision that we are going to consider in the following $+(z, x) = $-(z, x) 
and as a consequence the trapped surface has to be symmetric under a parity transformation 



3 We have to keep in mind that only D — 4 of the D — 3 coordinates d % parametrizing the sphere S D 4 are 
independent, since I'd 1 ) + . . . + (i? D ~ 3 ) 2 = 1. 
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in the longitudinal coordinate. This means that H + (r,9) = H-(r,9) = H(r,6), and we have 
that the symmetry conditions (3.10) imply that 



d e H(r,0)=d e H(r,7i)=0. 



(3.12) 



The numerical solution of the problem at hand is complicated by the fact that the function 
G{9) to be determined appears in the boundary conditions of the differential equation. Be- 
cause of this we would like to make a change of coordinates that transforms the complicated 
integration region into a simpler one. Following |12j . we define the new dimensionless radial 
coordinate p by 



r = LG(9)p. 



(3.13) 



Now, in terms of this new variable Eq. (|3.7|) reads 

2 



2 „2 



1 + G*p 



~G 



d >v^- 2 -p% d ^ + - 



(D - 3) + (D - 2)G 2 p 2 } d p 



P 

D-A 



~G 



G" 
~G 



D-A (G' 



tan# V G 



do — (D — 2)G 2 \ H(p 



(3.14) 



0. 



p 2 tan 9 

With this we are left with a boundary problem in the half-circle p = 1 with < 9 < tt, whereas 
the function G(9) defining the boundary in the original coordinates (r, 9) appears now in the 
equation. The boundary conditions on the lower and upper boundaries of the upper half-circle 
are respectively 



d e H(p, 0) = d e H(p, tt) = 0, d e H(p = 1,9) = — $ (r = LG(6), 9^j . (3.15) 
This means that we have von Neumann boundary conditions on the lower boundary and Dirich- 



let condition on the upper one. To derive the last equation in (3.15) we have used that the 
function defining the trapped surface vanishes on C, ^{p = 1,9) = 0, and as a consequence 
dg^/(p = 1,9) = 0. This equation actually simplifies the additional condition (3.11) 

G'(9) 2 



G(9Y 



AG(9Y 



(3.16) 



G(9) 2 

In order to implement the numerical algorithm it is convenient to introduce the dimensionless 
functions 9) = £tt(p,0), $(p,6) = z${p,0) and H(p,9) = \H(p,9). 
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Figure 3: Numerical solutions for the shape of the section of the surface C on the plane Y l -Y D 2 
in hyperbolic coordinates for D = 4,5,6,7 and 8 for various values of the impact parameter 



b and \i 
1 = 0-4 



b 



From the outer inwards the curves correspond to 



• 



' L 



0.2 



'' L 



0.5 (a 



and the critical value ^ 



The symbols indicate the position of the 



incoming sources. 



We can now attack the numerical solution of the boundary problem (3.14), (3.15) subjected 



to the additional condition (3.16). We follow the procedure used by the authors of [12] and 



choose a trial function G\{6) to solve the boundary problem (3.14), (3.15) using a method 



of finite differences. This gives a solution Hi(p,8) in terms of which we compute ^i(p, 0) = 
$(p, 6) — Hi(p,6). In the extremely unlikely case in which G\{9) were the solution of the 
problem the function 



U9) 



1 + 



9L 

G 



G 



0p*i(l, 



AG 2 



(3-17) 



would vanish. However, since this is not going to be the case we introduce a new trial function 

G 2 (9) = G 1 {e) + eT 1 {0), (3.18) 



with e > a small number, and solve again the boundary problem (3.14), (3.15) obtaining a 
second function H 2 (p, 6) and ^(p, #)• Computing now T 2 {6) as in Eq. (3.17) and iterating the 
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algorithm a sufficient number of times a numerical solution to the problem can be found with 
the required precision. 

We have used the procedure described above to find the solutions to the trapped surface 
equation in AdS^ with D = 4, 5, 6, 7 and 8. We have solved the boundary problem using a 
finite difference method with a grid of 50 points in the angular coordinate and 100 points in 
the radial one. With e = 10 -4 we find good convergence to the solution. In order to handle 
the coordinate singularity at r = we have imposed the von Neumann boundary condition 
d n H = at this point. 

Our numerical analysis shows the existence in all cases of a critical value of the impact 
parameter above which no solution to the marginally closed trapped surface equation is found. 
This is qualitatively similar to the corresponding situation in flat space [12J as well as in AdSs 
[TU] . In Fig. [3] we plot the section of the trapped surface on the Y l -Y D ~ 2 plane for a symmetric 
collision in D = 4, 5, 6, 7 and 8 with energy fi = G^L D ~ Z . For AdSs our numerical solutions 
are in perfect agreement with the results found by the authors of [10] using a different numerical 
method. 

On the left panel of Fig. |4j on the other hand, the values of the critical impact parameter 
as a function of the energy are shown for different dimensions. As can be seen in the right 
panel of Fig. |4| the results are very well fitted by the following scaling of the critical impact 
parameter with the energy of the incoming waves 



where the proportionality constant is of order one. One consequence of this scaling is that the 
dependence of the critical impact parameter with the energy, measured in units of the AdS/) 
energy L D ~ 3 G] V 1 , flattens as the dimension increases. In the limit D — > 00 (with L D ~ 3 G~^ 
fixed) this seems to implicate that the critical impact parameter is constant with the energy. 

In all our previous analysis we have focused our attention on physical situations where the 
energy of the incoming waves are larger or of the order of the AdS energy scale. This means 
that the physics is sensitive to the large-scale geometry of AdS. Flat space is retrieved by taking 
L — > 00 with the energy and impact parameter of the collision fixed. Since in this limit the 
energy is very small compared with the AdS natural energy scale, the flat space behavior of 
the critical impact parameter with the energy has to be recovered in the plots of Fig. [4] in 




(3.19) 
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Figure 4: In the left panel we plot the critical impact parameter ^ = b+ ~ L b ~ as a function of 9n-a for the 
symmetric collision of two shock waves with sources located at z = L. The curves correspond, from top to 
bottom, to D — 4, 5, 6, 7 and 8. On the right panel the same points are shown together with a fit to a function 



with the scaling of Eq. (3.191 



the region around the origin. On dimensional grounds, the scaling between the critical impact 
parameter and the energy for gravitational shock waves collision in flat space has to be of the 



form b c ~ (GV/z) D - 3 . A crossover between (3.19) and this scaling should take place as the 



energy of the collision is much smaller than the AdS energy scale. 

To compare our results with the analytical study of the formation of trapped surfaces in 
the collision of shock waves in AdSs with nonvanishing impact parameter carried out in [8], 
we notice that the authors of this reference work in the limit where the energy of the collision 
satisfies 



L D - 3 



D-2 



> 1, (3.20) 



while keeping fixed the size of the impact parameter. Because of the scaling (3.19) the results 
of [8] are in the regime where the impact parameter is always much smaller than the critical 
value and, as a consequence, a marginally closed trapped surface always forms. 



4 Discussion 

Having obtained the numerical solution for the collision of two shock waves with "spatial" 
impact parameter we can apply the transformations described in the previous section to get 
solutions of the closed trapped surface equation for the collision of two waves with "holographic" 
impact parameter. On the left panel of Fig. [5] we have plotted, for various values of the spatial 
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JC'/L 




x/L 



Figure 5: In the left panel we plot the section of the closed trapped surfaces in Poincare coordinates (x' , z') 



shown in Eqs. (2.9) and (2.26). The curves, from outer inwards, correspond to the following values of the 
impact parameter: \ = (•), \ = 0.3 (■), \ = 0.6 (♦), \ = 0.8 (A) and the critical value % (▼). On the 
right panel we see the sections of the closed trapped surfaces after inverting the rotations (2.9) back to the 
coordinates (x,z). In both plots the horizontal axis representes the boundary of H3. 



impact parameter, the section of the marginally trapped surface for a symmetric collision in 
ACIS5 of two waves with energies n± = GnL~ 2 . Unlike the plots shown in Fig. [3] here we have 



used the Poincare coordinates (x',z') defined in Eq. (2.20). Hence, in this plot the boundary 



of AdSs is located along the horizontal axis. On the right panel of Fig. |5j on the other hand, 
we see the section of the trapped surfaces previous to the 0(2) rotation of angle | described 
in Section [2j Now the sources lie in the z axis and therefore the contours show the section 
of the trapped surface formed as a result of a collision of two waves with holographic impact 
parameter. The innermost curve represents the surface corresponding to the critical value of 
the impact parameter, above which no marginally closed trapped surface of the type studied 
here exits. 

The conclusion is the existence of a maximum value of the holographic impact parameter 
for the formation of this class of trapped surfaces. Taking this fact at face value, one could be 
tempted to interpret it holographically as signaling a threshold for the formation of a thermal 
CFT plasma. Thus, thermalization would not take place when the size difference between the 
two colliding energy lumps is larger than the critical value. This would actually make sense 
on physical grounds: if one of the colliding lumps is much smaller than the other it might not 
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have enough degrees of freedom to induce a thermal state in the whole system as a result of the 
collision. One should be reminded, however, that the non existence of a closed trapped surface 
of the type studied above does not rule out the possibility of other trapped surfaces being 
formed in the process. Therefore, in order to decide whether a threshold for thermalization 
exists as a function of the difference in size of the colliding objects one would need to solve for 
the geometry of the interaction region of Fig. [TJ This problem, together with other issues such 
as the thermalization time for subcritical collisions, will be addressed elsewhere. 

In the dual conformal field theory, the 0(2) transformation that maps a collision of two 
equal objects with impact parameter into the head-on collision of two lumps of different size is 
a combination of a translation and a special conformal transformation in one of the transverse 



coordinates, inherited from the AdS.o bulk transformation (2.10). Interestingly, this relates 
a physical process with nonvanishing angular momentum to another one where the angular 
momentum is zero. In the gravitational dual, on the other hand, the combination [see Eqs. 



(2.14) and (2.17) 



Q± = /4(0)£ 



1 + 



z' ± (ey + b' ± (e) 



(4.1) 



L 2 

is independent of the angle 9 and therefore invariant under the 0(2) rotation. Holographically, 
this relates the impact parameter, the size of the lumps and their energies. A second invariant 
can be also constructed as 

Q{9) = Li> + W_{e)z' + {e)z'_{9) (4.2) 

that mixes the energies of the colliding objects with their sizes. 

In D = 5 the collision of two gauge theory energy lumps with large size difference could be 
used as a model for hadron-nucleus collisions at strong coupling. Here we have considered the 
case of head-on scattering. Collisions of unequal objects with a nonvanishing impact parameters 
can be studied as well by performing 0(2) rotations on the numerical solutions obtained in the 



previous section with an angle \9\ < §. It would be interesting to see if the invariants (4.1) and 



(4.2) have any relevance in the phenomenological description of this type of collisions. 
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Appendix 



In this Appendix we give some details on the change of coordinates that have to be used to 
eliminate the delta-function terms in the metric of the shock wave and therefore to make the 
congruence of null geodesies continuous across the wave fronts u — 0, v — 0. Let us begin with 



the metric (1.1) now written in the form 

L 2 



dz — dudv + dx + <p+(z, x)5(u)du 



(A.3) 



and introduce the new coordinates (U, V, Z, X) defined by 
U. 



u 



x 



V + <f + (Z, X)6(U) + \U6{U) [v<p + (Z, X)] 2 + \U6{U) [d z <p+(Z, A)] 
Z+±W{U)dz<p+(Z,X), 

x + lue(u)v<p+(z,x). 



(A.4) 



The presence of the step function 6(U) in the change of coordinates is necessary in order 
to eliminate the delta-function terms in the metric, responsible for the discontinuity of the 



geodesies. After a long but simple calculation we find that the metric (A.3) for the first wave 
takes the form 



ds< 



-dUdV + 'H { X } 'H { ^ ) dX a dX c 



z + wo(u)d z <p A 



-i 2 



(A.5) 



with 



H 



(+) 

ah 



(I 



1 



ab 



U9(U)d a d b( p+(Z,X). 



(A.6) 



The indices a, b, c run over all the spatial coordinates (Z,X). The line element (A.5) is the 
AdS analog of the so-called Rosen form of a plane wave in flat space-time, whereas the metric 



(A.3) we started with is the Brinkmann form of the wave. 
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The analysis can be repeated for the second wave, whose metric can be obtained from 



Eq. (A. 3) by changing u — > v and (f + (z,x) — > (p~(z,x). Implementing then the change of 



coordinates (A. 4) with the replacements U — > V, u — > v and <p+(Z, X) — > (p_(Z,X) we arrive 
at the metric 



dUdV + ■H < a ~ ) n i ~ ) dX a dX 



z + \ve(y)d z y 2 



(A.7) 



where now the metric function %L is given by 



U 



(-) 



5 ab + -V6{V)d a d m {Z,X). 



(A.i 



ab -uu . 2 

It is important to notice that the metric elements ( |A.5 ) and ( |A.7 ) coincide in the lower wedge 
U < 0, V < where the line element of AdS^ is recovered. Hence, the metric in the whole 
region outside the interaction wedge U > 0, V > can be written as 



-dUdV + 



ds 2 



nj(+)nj(+) I nj(-)nj(-) 

^ab n-bc n-ab ^bc 



ab 



dX a dX c 



(A.9) 



z + \u6{u)dzv + + lve(y)d z <p- 

The next step is to write the Einstein equations in the new coordinates. This task, however, 



is simplified by noticing that for the line element of the first wave (A. 3) the only nonvanishing 
component of both + Ag^ and is the u-u component. Since U = u we find again that 
the equation for <p+(Z, X) is given by Guu + Agjju = StiG^Tjju and after a bit of algebra Eq. 



(2.4) is recovered in the new coordinates (the same argument applies to the second wave with 
the replacement U — Y V) 

D-l 



□ 



D-2 



H D _ 2 



$±(Z,X) = -167rGV/i ± 



5(Z -z±)5 (D - 3 \X -b±), (A.10) 



where $±(Z,X) = ^(p±(Z,X). In fact, the reason behind obtaining the same equations in 
both coordinate systems is that the change (z,x) — > (Z,X) is trivial in the transverse space 
M.d-2- This can seen inmediately by setting U = in (A. 4), and respectively V = in the 
corresponding change of coordinates for the second wave. 
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